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The X-ray nebula of the filled center supernova remnant 3C58 
and its interaction with the environment"^ 



o 
o 

(N 

D 
tin 



> 
00 

o 

O 

o 

:^ 

o 
U 



F. Bocchino^'^, R. S. Warwick^, P. Marty"^, D. Lumb^ W. Becker^ and C. Pigot*^ 

Astrophysics Division, Space Science Dept. of ESA, ESTEC, Postbus 299, 2200AG Noordwijk, The Netherlands 
Osservatorio Astronomico di Palermo, Piazza del Parlamento 1, 90134 Palermo, Italy 
Department of Physics and Astronomy, University of Leicester, Leicester LEI 7RH, U.K. 
Institut d'Astrophysique Spatiale, Campus Universite Paris-Sud, F-91405 Orsay Cedex, France 
Max-Planck-Institut fiir extraterrestrische Physik, Giessenbachstarsse 1, D-85740 Garching, Germany 
DAPNIA / Service d'Astrophysique, CEA/Saclay F-91191 Gif-Sur-Yvette, CEDEX, France 

Received 5 December 2000; Accepted: 30 January 2001 

Abstract. An XMM-Newton observation of the plerionic supernova remnant 3C58 has allowed us to study the 
X-ray nebula with unprecedented detail. A spatially resolved spectral analysis with a resolution of 8" has yielded 
a precise determination of the relation between the spectral index and the distance from the center. We do 
not see any evidence for bright thermal emission from the central core. In contrast with previous ASCA and 
Einstein results, we derive an upper limit to the black-body 0.5-10 keV luminosity and emitting area of 1.8 x 10^^ 
erg s~^and 1.3 x 10^" cm^, respectively, ruling out emission from the hot surface of the putative neutron star 
and also excluding the "outer-gap" model for hot polar caps. We have performed for the first time a spectral 
analysis of the outer regions of the X-ray nebula, where most of the emission is still non-thermal, but where the 
addition of a soft (kT=0.2-0.3 keV) optically thin plasma component is required to fit the spectrum at i? < 1 
keV. This component provides 6% of the whole remnant observed flux in the 0.5-10.0 keV band. We show that 
a Sedov interpretation is incompatible with the SN1181-3C58 association, unless there is a strong deviation from 
electron-ion energy equipartition, and that an origin of this thermal emission in terms of the expansion of the 
nebula into the ejecta core nicely fits all the radio and X-ray observations. 
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1. Introduction 

3C58 is a beautiful example of a fiUed-center (or pleri- 
onic) supernova remnant (SNR), probably associated with 
the historical supernova event in 1181 A.D. ( Stcphensoii| 



1971). lihis object has always received much attention be- 
cause, on the one hand, it shows some characteristics sim- 
ilar to the Crab SNR, while, on the other hand, it seems 
very different to the Crab itself. For instance, 3C58 has a 



point source located in the center, as would be expected, 
since the morphology strongly suggest that the nebula is 
powered by a spinning neutron star. Despite considerable 
effort, pulsations have not been detected so far in either 
radio (see e.g. Lorimer et al. 199S), nor in X-rays (Helfand 



et al. 1995, H95 hereafter). Moreover, the X-ray to radio 



compact (iu' X p') elliptical morphology with a very bright 



flux ratio (fx/fr) of 3C58 is 100 times lower than that of 
the Crab (H95), its spectral break occurs at 50 GHz (300 
times less than the break of the Crab, [Green & Scheuei 



core (Reynolds k, AUer iy»»j and linear size similar to the 
Crab; Frail & Moffett (1993) have reported a wisp-like 



1992) and its radio luminosity is increasing instead of de- 
creasing as expected ( preen 1987 ). These remarkable dif 



elonga jed structure at z.p ' from the core, which has been 
observed also m the Crab and which is probably associ- 
ated with the pulsar wind termination shock. However, 
unlike the Crab, there is no clear evidence of a pulsating 



ferences are also seen in other plerions, and Woltjer et al 
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(1997) have proposed the introduction of a new sub-class 
of plerions, the "non Crab- like plerions", of which 3C58 
can be considered the prototype. For these objects, a non- 
standard evolution of the pulsar can be invoked, but the 
details are not yet clear. 

It is therefore very important to investigate the phys- 
ical properties which render 3C58 so peculiar, in order to 
put this object and its sub-class in the right perspective. 
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In particular, the detection or non-detection of the cen- 



tral source is obviously a key point. Becker et al. (1982) 
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3C58 source profile 



reported the presence of a compact X-ray source in 3C58 
from their Einstein HRI observation, about 10" in extent, 
and contributing to 15% of the detected X-ray flux. Later, 
Helfand et al. (1995) revisited the X-ray emission of 3C58 
using ROSAT HRI data, confirming the compact source 
and favored a model in terms of hot polar caps to explain 
the emission. However, it has not been possible so far to 
take an X-ray spectrum of the source to study it, and to 
understand if it is rea lly a point source o r an enhancement 
of the pulsar nebula. Torii et al. (2000) have pointed out 
that the inclusion of a black-body component in the fit of 
the ASCA GIS-I-SIS data of 3C58 yields an improvement 
of the x^. They claim that the best-fit black-body com- 
ponent is responsible for ~ 7% of the unabsorbed flux in 
the 0.5-10 keV of the whole remnant, and that it is the 
spectral signature of the central source. 

It is also very important to assess the presence of a 
shell around the pulsar nebula, for it may give compelling 
constraints on the age of the remnant, the shock veloc- 
ity and the density of the environment. In the case of 
the ~ 800 yr old 3C58 (as other plerions as well), it is 
expected that the main shock will encounter the stellar 
ejecta and/or the interstellar medium (ISM) giving rise 
to a limb brightened shell. However, no sign of a shell 
has been detected at centimeter wavelengths at distance s 
greater than 5' from the core ([Reynolds fc AUer 1985| ). 



However, Reynolds fc AUer (1988) have imaged the faint 



outer emission of the nebula at a distance between 2' and 
4' and have noticed limb brightening at several locations. 
In this paper, we present a study of the XMM-Newton data 
of 3C58 obtained during the Calibration and Performance 
Verification (Cal/PV) phase of the mission. We show that 
a stringent upper-limit can be imposed on the presence 
of thermal emission from a central point source, and we 
also find evidence for a thermal outer shell. The implica- 
tions of these findings for current models of 3C58 are also 
discussed. 



2. Data analysis 

2.1. Observations 

3C58 has been observed as part of the Cal/PV phase of 
the XMM-Newton satellite ( [lansen et al. 20001 ). In this 
paper, we focus on the 25 ks on-axis observation performed 
during orbit 47, on March 12th 2000. Data from the two 



MOS (Turner et al. 2000) cameras and the PN (Strueder 



et al. 200C) camera have been used. MOS and PN cameras 
are CCD arrays which collect X-ray photons between 0.1 
and 15 keV and have a field of view of 30'. The pixel size 
is 1.1" and 4.1" for MOS and PN respectively, and this 
should be compared with the mirror PSF width of 6" — 
15" FWHM-HEW. The data have been acquired with the 
medium filter and in full image mode, and therefore the 
temporal resolution is low, 2.5 s and 73 ms for MOS and 
PN, respectively. The worse spatial resolution of the PN 
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Fig. 1. MOSl+2 source profile, centered at 2''05"38.0'' 
-|-64''49™37'' (J2000V Soft, medium and hard energy band 
as defined in Fig. g are represented in red, green and 
blue, respectively. The solid black line is the Point Spread 
Function (PSF) of the mirror module 3 at 1.5 keV 



(Aschenbach et al. 2000) 



is compensated by its greater sensitivity, on the average 
20-30% more than the combined two MOS. 

The Standard Analysis System (SAS) software we have 
used (version 5.0 alpha, xmmsas-20001011-1559) takes 
care of most of the required events screening. However, we 
have further screened the data to eliminate some residual 
hot pixels and occasional background enhancement due 
to intense flux of soft protons in the magnetosphere. In 
particular, we have extracted the background light-curve 
from a region free of sources, and we have identified time 
intervals of unusually high count rates (more than twice 
the "quiescent" background rate) and removed them from 
subsequent analysis. Moreover we have also selected events 
with value of the PATTERN pipeline assigned keywords 
between and 12. The exposure time of the screened ob- 
servations is 16 ks for the MOS and 12 ks for the PN 
detectors. 

2.2. X-ray morphology 

In Fig. 1^, we report the 3C58 profile as seen by the 
MOSl-l-2 computed in 65 annuli with Ar = 4" up to a 
distance of ^ 4' from the centroid, located at 2''05'"38.0'' 
-|-64''49™37'' (J2000, ±1"). The background was collected 
from a ring with Ar = 30" located at 4.5' from the center, 
and has been chosen to fall entirely in the central MOS 
chip. The comparison with the expected mirror PSF, also 
shown in Fig. ^ indicates that the central source is ex- 
tended. The bright core region can be defined up to ~ 50" 
from the center; Further out, the slope changes and the 
behavior of the soft curve and hard curve is different. The 



X-rays from 3C58 
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hard curve declines with the same slope down to the limit 
of the radio nebula at ^ 200", while the soft curve is 
clearly flatter than the hard one, with a sudden change of 
slope occurring at 140". 

The core is elongated in the N-S direction, with a 
FWHM of 13.5 ±1.5" in this direction (9.0±1.5" in the E- 
W direction), in agreement with the ROSAT HRI results 
presented by H95. The centroid of the MOSl emission is 
located 4" and 8" south of the ROSAT and Einstein po- 
sitions. The data confir m that the radio wisp observed by 
Frail fc Moffett (1993)| is located to the East of the X-ray 
peak, at 5" from the the MOSl centroid. However, the 
uncertainties in real attitude reconstruction, which is still 
not fully implemented, may give an uncertainty in abso- 
lute position determination as large as ~ 5". 

Fig. H shows the summed MOSl+2 images of 3C58 in 
three different bands, namely 0.1-1 keV, 1-2 keV and 2-10 
keV bands, along with the weakest 1446 MHz radio con- 
tour at 0.15 Jy beam"'^ and the contour at 0.6 Jy beam~^ 
of Fig. 3 in [Reynolds fc AUer (1988)| (both in white in 
Fig. H). A small separation of the two contours (e.g. in 
the SW) indicates limb brightening and therefore confine- 
ment i n the source. In addition, we also shown the X-ray 



8 annuli, we have also defined a region (the "edge" re- 
gion hereafter), represented by the union of two ellipses 
matching the outer X-ray edge of the nebula|^, minus a 
circle with same radius as the 8th annulus. This region 
is particularly suited for the study of the X-ray emission 
coming from the outermost fringes of the X-ray and radio 
nebula. 

We have used three different emission models to fit the 
3C58 data, namely a power-law model, a power-law model 



plus the optically thin plasma model of Mewe et al. (1985) 
with Fe L calculation of Licdahl et al. (1995), pl-|-MEKAL 
hereafter, and a power-law model plus a black-body spec- 
trum, pl-|-bbody hereafter. These three models encompass 
what we could possibly expect from an X-ray nebula, the 
last two representing eventual contributions from a ther- 
mal shell (as in the case of known plerion-composite SNR), 
and from a compact source in the center as pointed out by 



Helfand et al. (1995). The temperature of the pl-|-MEKAL 



model and of the pl+bbody models have been constrained 
in the 0.1-10 keV and i n the 0.1-2.0 keV, respectiv ely. 
Abundances are those of Anders & Grevesse (1989). All 



contou 's (in yellow), and it is clear that there is corrcspon- (1983) 



the models have been modified b y the interstellar absorp- 
tion ac cording to cross-sections of Morrison fc McCammon 
where we have let the equivalent hydrogen col- 



dence between the radio contours and the X-ray emission 
of the soft image. In particular, the weakest radio contour 
matches the outer edge of the weak extended X-ray emis- 
sion, especially at the north and south edges. Fig. |^ also 
shows that the size of the nebula de creases as the ener gy 
increases, an effect also reported by Torii et al. (2000)| . 



umn density Nh vary. Since we have noted that the non- 
thermal component provides most of the flux in the XMM- 
Newton bandwidth, and that the residual thermal compo- 
nents of the model pl-|-MEKAL and pl-t-bbody are corre- 
lated with the value of Nh , we have also performed a set of 
fittings fixing the Nh value to 4x 10^^ cm'^, which is com- 



patible with previous estimates (Helfand ct al. 1995, Torii 



2.3. SpecLidl diidlybib 

The high spatial and spectral resolution of the instru- 
mentation aboard the XMM-Newton satellite allows us 
for the first time to perform spatially resolved spectral 
analysis of the X-ray nebula associated to 3C58. To this 
end, we extracted spectra from 8 concentric annuli with 
Ar = 8" centered at the same position of the X-ray 
profiles of Fig. |l|, covering the core emission up to 1.1' 
from the center. Given the uncertainties in present cal- 
ibration of the MOS and PN cameras, and the energy 
dependence of the vignetting correction above 5 keV, we 
have restricted our spectral analysis to photons in the 
0.5-5 keV energy band. The spectra have been back- 
ground subtracted using the same background region in- 
troduced in the previous subsection, and have been re- 
binned to ensure that a minimum of 30 counts are present 
in each energy channel. As for response matrices and 
effective area files, we have used the latest version of 
standard MOS and PN matrices provided by the cal- 
ibration team (mosl_medium_all_qel7_rmf3_tel5_15 and 
epn_fs20^Y9 -medium). We have summed the spectra ex- 
tracted with MOSl and M0S2, and we have also rescaled 
the response matrix to reflect this operation^. Besides the 



et al. 2000) and it is also consistent with the result we 



obtained letting it vary. 

2.3.1. The 7 vs. radius relation 

Fig. H shows the best- fit value of the absorption and of the 
power-law photon index (7) as a function of the distance 
from remnant center, obtained with a fit to a power-law 
emission model only. The data clearly show the effect of 
synchrotron burn-off of high energy electrons as the radius 
increases. This effect has also been observed in G21.5-0.9 
both with Chandra ([Slane et al. 2000D and XMM-Newton 
(Warwick ct al. 2000) and it is related to inhomogene- 
ity in the particle distribution inside the plerion nebula. 
The straight line in the lower panel of Fig. || represents 
the linear best-fit to the 7 — r relation, j = A + Br with 
A = 1.97 ± 0.03 and B = 8.0± 0.6 x 10"^ and r is in arc- 
seconds. The MOS fits to the single power-law reported in 
Fig. H are statistically acceptable^ from ring 1 to ring 6, 
while not acceptable in rings 7-8 and in "the edge" region. 
This is true for both Nh free and Nh fixed fits, and Table 



^ We will refer hereafter to the summed spectra as simply 
the MOS spectra. 



The E-W and N-S ellipses have large and small semi-axis 
of 3.3' X 1.1' and 1.7' x 1.1', respectively. 

^ We regard a fit as statistically acceptable when the null hy- 
pothesis probability computed from the reduced a-nd num- 
ber of degrees of freedom is greater then 5%. 
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Table 1. Goodness of MOS fits to the 8 annuli of 3C58. The absorption is fixed to 4 x 10^^ cm^^. For each model, we 
report the value of x^, degrees of freedom (dof), null hypothesis probability. For the pl+MEKAL and pl+bbody model, 
we also report the value of the best-fit temperature. The last column is the unabsorbed flux in the 0.5-10.0 keV in 
units of 10~^^ erg cm^^ s^^; its statistical error is ±0.1. 



Reg. Power-law Pow-|-mekal Pow-|-bbody 

X^/dof Prob. kT X^/dof Prob. kT X^/dof Prob. flux 



1 


74/73 


46 


1.0(> 0.4) 


71/71 


48 


0.2(> 0.1) 


73/71 


40 


15.3 


2 


69/65 


35 


6.0(> 0.1) 


67/63 


35 


0.9(> 0.1) 


68/63 


30 


13.0 


3 


65/63 


41 


0.6(> 0.1) 


65/61 


34 


0.5 ±0.3 


62/61 


44 


11.4 


4 


52/59 


71 


0.1(> 0.1) 


51/57 


70 


0.1(> 0.1) 


51/57 


69 


12.3 


5 


54/61 


72 


0.7(> 0.1) 


52/59 


72 


0.2(> 0.1) 


53/59 


69 


12.0 


6 


52/59 


74 


0.1(> 0.1) 


50/57 


72 


0.6(> 0.1) 


50/57 


72 


10.9 


7 


96/60 


0.2 


0.2(< 0.25) 


91/58 


0.4 


0.81°:? 


91/58 


0.4 


10.7 


8 


79/58 


3 


'-'■''-0.1 


71/56 


8 


0.81°;? 


73/56 


7 


10.6 



|l| reports the x^/do/ values of the fits. It is interesting to 
note that the Njj value, when left free to vary, is signifi- 
cantly lower than the average value of 4 x 10^^ cm"^ for 
the outer nebula regions (Fig. ^). Moreover, the best-fit 7 
of the shell is off the trend dictated by fit to the spectra 
of the rings when Nh is left free to vary, while it shows 
lower deviation in the fit with Nh fixed. If the X-ray emis- 
sion of the outer rings is dominated by the non-thermal 
component of the plerion, we do not expect significative 
variation of the absorption, and the data seem to confirm 
that a fixed Nh may be more appropriate. However, the 
fact that fits of the outer rings and "edge" data are less 
acceptable than fits to inner rings data strongly suggests 
a contribution from other components. 

2.3.2. Additional thermal components 

In order to assess the presence of any additional emission 
from a central compact source (as suggested by II95) and 
from any thermal shell, we now consider the results of fits 
with pl±MEKAL and pl±bbody models. The inclusion of 
the additional component (mekal or bbody) does not in- 
crease the null hypothesis probability above 5%, except for 
ring 8. This is also reported in Table |l|, in which the reader 
finds the data needed to evaluate the goodness of the fits 
we have performed. However, it should be noted that the 
inclusion of the mekal component leads to a reduction 
of x^ (and therefore to a better fit in a relative way) for 
the outermost rings, and a dramatic reduction in case of 
the "edge" region, specially for MOS. This is also shown 
in Fig. 1^, where we show the MOS "edge" spectrum along 
with its power-law only and pl±MEKAL best-fit model and 
residuals, and in Table |^ where we summarize the results 
of spectral fitting to the "edge" spectrum. In Fig. ||, we re- 
port the ratio of the fiux of the second component (/mckai 
from the pl±MEKAL fit in the upper panel, and /bbody 
from the pl±bbody fit in the lower panel) to the total fiux 
of ring 1-8 and the "edge" region. 

The fit with an additional black-body component to 
ring 1 places an upper limit to the unabsorbed flux due to 



this component of 20% of the total 0.5-2.0 keV flux in this 
region (i.e. 1.5 x 10""'^'^ erg cm~^ s~^, or a luminosity Lx = 
1.8 X lO'^^I?! 2 srg s~^, where 1)3.2 is the distance in units 
of 3.2 kpc, which i s the most reliable value according to 
Roberts et al. 1993[) and 10% of the total 0.5-10.0 keV flux 



of this region (i.e. 1.6 x 10"^'^ erg cm~^ s~^, same Lx as 
before). Since the whole remnant 0.5-10 keV flux is 1.6 x 
10~^^ erg cm~^ s~^, the upper limit corresponds to 1% 
of the remnant fiux, significantly below the 7% found by 
Torii et al. (2000). Fig. |^ also shows that spectral fittings 
of rings 2-6 yields only upper-limits to the presence of an 
additional mekal or black-body thermal emission. 

On the other hand, the situation is different for ring 8 
and the "edge" region (Table |]and||), where the addition 
of the thermal component yields a significant decrement 
of the x^ (according to an F-test) and the thermal flux 
contribution to the 0.5-2.0 total flux is between 5% and 
30%. In particular, the pl±MEKAL fit in the "edge" region 
suggests that between 10 and 20% of the 0.5-2.0 keV flux 
from this region (5—9 x lO^^'^ erg cm~^ s~^, or Lx = 5.8 — 
10.0 X lO'^^I?! 2 ^^S is due to the thermal component 
(8-14% if we consider the 0.5-10 keV band). If compared 
to the whole remnant 0.5-10 keV flux, the thermal soft 
excess yields a contribution of 3%-6%. It is important to 
note that in this region the pl±MEKAL model is to be 
preferred over the pl±bbody model (Table ||). 



2.3.3. PSF effects 

We have also investigated whether our results are affected 
by any effect related to the Point Spread Function (PSF) 
of the XMM-Newton mirror. Since only 55% of the pho- 
tons of a point source are contained within a circle of 8" 
radius, we expect a certain amount of scattering outside 
the rings. However, this effect is not likely to modify our 
conclusions. In the central ring, for instance, the strongly 
peaked morphology keeps the contamination from ring 2 
low. Moreover, the scattering can only smooth the 7 — r 
relation in Fig. ||, so the observed slope is, strictly speak- 
ing, a lower limit to the real slope. Finally, in order to 
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Fig. 4. A close up of MOS spectrum for the "edge" region. Left panel: the fit with a simple power-law model does not 
well represent the spectrum below 1 keV. Center panel: the pl-|-MEKAL model with kT~0.2 keV improves the description 
of the data. The spectral region near the oxygen edge (between 0.4 and 0.6) is know to have some residual calibration 
problem. Right panel: Unfolded best-fit pl+MEKAL model. Dot-dashed line represents the thermal component, which 
shows emission lines and falls off rapidly above 1 keV, where the power-law component dominates. 



Table 2. PN and MOS fit results of the "edge" region of 2.4. Timing analysis 



3C58. All the fits yield an absorbed flux of 6.1 x IQ-^^ erg 



cm 



in the 0.5-10 keV band. 



Model 



kT 
keV 



xVdof 







MOSl±2 






pow-law 


3.5 ±0.2 


2.61 ±0.06 




245/170 


pl-l-MEKAL 


7 O+0.6 
' •■'-0.9 


2.89 ±0.08 


0.18 ±0.01 


202/168 


pl-|-bbody 


i; c+0.5 
^•^-0.8 


2.75 ± 0.06 


< 0.11 


231/168 


pow-law 


4.0 


2.75 ± 0.04 




261/171 


pl-l-MEKAL 


4.0 


2.61 ±0.05 


0.25 ± 0.03 


226/169 


pl-|-bbody 


4.0 


3.07 ±0.07 


0.78 ± 0.08 


243/169 






PN 






pow-law 


3.1 ±0.3 


2.42 ±0.11 




113/138 


pl-|-MEKAL 




9 71-1-0.15 

^■'-■■-o.ii 


0.18 ±0.03 


107/136 


pl-l-bbody 


5.2tl:? 




< 0.13 


109/136 


pow-law 


4.0 


2.66 ± 0.06 




130/139 


pH-MEKAL 


4.0 


2.47 ±0.10 


0.24 ± 0.06 


110/137 


pl-hbbody 


4.0 


2.46 ±0.10 


< 0.13 


110/136 



verify whether the soft component detected at large radii 
is an artifact of the PSF. wc; have followed this procedure: 
we have fitted the summed spectrum of rings 1-8 with a 
power-law model, and we have found Nh = 3.9 x 10^^ 
cm~^ and 7 ~ 2.21; then, we have fitted the spectrum of 
the "edge" region with the best-fit model found for ring 1- 
8 with fixed Nh and 7, plus a second absorbed power-law 
model with free Nh and 7. The normalizations of the two 
models were left free to vary. We have found that the best- 
fit normalization of the "fixed" power-law is zero, suggest- 
ing that negligible contamination from the spectrum of the 
ring 1-8 is observed in the "edge" spectrum. 



In order to search pulsed X-ray emission from the cen- 
tral source which is powering 3C58, we have calculated 
the power spectrum density of the time series of events of 
the central ring used for spectral analysis. Unfortunately, 
the time resolution of the EPIC cameras when operated 
in full image mode is low, and we were not able to sam- 
ple frequencies above 1 Hz. For timing analysis, wo have 
collected MOSl and PN photons in ring 1 between 0.5 
and 10 keV, and we have analyzed the two independently. 
Furthermore, we have also analyzed a more restricted en- 
ergy range, 2-10 keV. 

None of the power spectra show features significant at 
the 99% confidence level. The corresponding MOSl upper 
limits to relative amplitude of a sinusoidal pulsed signal 
in the 5 x lO'^ - 0.1 Hz is 6.2% and 10.0% for the 0.5-10 
keV and 2-10 keV respectively, while the PN upper limits 
in the lO'^ - 1 Hz are 2.0% and 3.1% in the broad and 
hard band, respectively. 

3. Discussion 

The good spatial and spectral resolution XMM-Newton 
data on 3C58 has allowed us to perform a detailed spa- 
tially resolved spectral analysis of this remnant, and in 
particular, to address the topic of the spectra of the com- 
pact core at its center and the presence of an X-ray shell. 

3.1. The compact core 

The possible emission mechanisms for the compact core 
of 3C58 have been reviewed by H95, and a preference has 
been given to a thermal model for the emission from hot 
polar caps, mostly on the basis of a review of ROSAT 
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and Einstein HRI data. Torii et al. (2000) have confirmed 
the detection of an additional component in the remnant 
spectrum, which they have modeled with a black-body of 
kT=0.40 - 0.45 keV and a flux of 1. 



erg s 



in the 0.5-10 keV band. While the [Torii et al. (2000)| flux 
estimate is within 20% of the H95 one, we have found in 
ring 1 an upper limit to the flux of the kT=0.25 keV black- 
body component a factor of 10 less than the H95 estimate. 
The inferred upper limit on the emitting area is 1.3 x 10^" 
cm^ (corresponding to a radius of ~ 650 m), more or less 



the value found by Torii et al. (2000) , and this is expected 



since our best-fit temperature is a factor of two less than 
the one deri ved by ASCA. By f orcing the temperature at 
the value of [Torii et al. (2000)| , the corresponding upper 
limits are 9.4 x 10"" erg cm^^ g-i ^ ^.1 x lO^'^Dl ^ 
erg s~^) and 8 x 10* cm^. While this large difference in 
the inferred black-body flux may seem puzzling, it may 
be explained by the dramatic improvement of the XMM- 
Newton spatial resolution compared with the ASCA one. 
In fact, if we compute the sum of our best-flt values for 
the flux of the black-body component of the pl+bbody fit 
over all the rings and the "edge" region, we obtain a 0.5- 
10 keV flux of 1.2 x 10~^^ erg cm~^ s~^, or a luminosity 



of 1.3 



erg s ^, in very good agreement with 

More 



the value reported by H95 and Torii et al. (2000) 



than 50% of this flux comes from ring 7-8 and the "edge" 
region: it is therefore probable that the ASCA data cor- 
rectly detected a soft excess in the integrated spectrum of 
3C58, but this was incorrectly attributed to the thermal 
component of the central compact object, whereas Fig. 
^ suggests that it comes from the external parts of the 
nebula. The more stringent XMM-Newton upper limit to 
the black-body component also heavily constrains its in- 
terpretation. Our value of the emission area is too low to 
assume that the whole surface of neutron star is loosing 
the residual heat of formation, and therefore would sug- 
gest the hot polar caps mechanisms. A review of the polar 



caps heating mechanisms can be found by Yancopoulos 
et al. 1 1994) . H95 favors the "outer-gap" model, and in 



this case we expect Lx = IQ Bi2P^ erg s~ , where B12 
is the magnetic field in units of 10^^ G and P is the period 
in seconds. As shown in Fig. ^, our Lx upper limit yields 
a magnetic field higher than the highest value known for 
allowed B and P values of an 810 yr old pulsar. For this 
model to be compatible with the Lx upper limit, either 
an older pulsar must be present or, unlike the Crab, the 
pulsar has not suffered high energy losses, in which case 
all the points below the isochrone lines of Fig. ^ are al- 
lowed. Other mechanisms invoked to explain the hot po- 
lar caps yield estimates of Lx lower then the "outer-gap" 
model and in principle may be compatible with the obser- 
vations. However, other observed pulsar nebulae-thermal 



pulsar pairs reported by [Seward fc Wang (1988) show a 
ratio between compact object luminosity and total (pul- 
sar -|- nebula) luminosity Lc/itot of ~ 0.1 — 0.2 (the Crab 
has Lc/itot = 0.04), while, using the observed total 3C58 
luminosity of 1.8 x lO^'' erg s~^and the upper limit we 
have derived, we derive a ratio < 0.01. Fig. ^ also shows 



that we expect a period less then 1 s for reasonable val- 
ues of B, and therefore, to properly detect a pulsation, 
we would need a better time resolution than those of the 
EPIC camera in full image mode. 

3.2. Thermal emission at the edge of the nebula 

The interpretation of the additional component in terms 
of expansion of the main shock in the environment of the 
SNR cannot be excluded a priori: recently, there have 
been some cases of SNRs recla ssified as com posite (e.g . 
Gil. 2-0.3, [Vasisht et al. 1996| ; G327.1-1.1, [Sun et al 



1999 ). The observed temperature of the mekal compo- 



nent corresponds to a shock speed of ~ 450 km s , 

while the emissivity corresponds to a post-shock density 

of (0.60 ± 0.15)1)3 cmT^ and a pre-shock density four 

times smaller (assuming the emission comes from a thin 

shell at r = 2.5', where the soft profile in Fig. |l| drops 

rapidly, and Ar — r/12). The emitting mass is of the or- 
5/2 

der of O.ID3 2 M0. The inferred density is very reasonable 
for typical low galactic latitude Sedov SNRs, but on the 
other hand it can hardly be reconciled with the associa- 
tion between SN1181 and 3C58. This is shown in Fig. 0, 
which reports the allowed values of the SNR age and of 
the distance according to a simple Sedov analysis following 
the outline given in Kassim et al. (1993). It is clear that 
from a pure geometrical point of view, the measured X- 
ray temperature may represent a Sedov shock of an ~ 800 
yr old SNR located at a distance of 3.2 kpc ("Geom." 
solid line in Fig. |^), but the measured emission measure 
of the MEKAL component yields a solution with a distance 
well above 10 kpc and and a remnant age well in excess 
of ^ 10** yr, if we assume that the explosion energy is of 
the order of the "canonical" value 10^^ erg (logE'si — 
in Fig. ^, where £'51 is in units of 10^^ erg). Among the 
lowest values of the explosion energy quoted in the liter- 



ature, we have £^51 = 0.1 for Vela (Bocchino et al. 1999) 
and 0.02-0.3 for G292.0-M.8 ( [Hughes fc Singh 1994[ ), arid 
we note that, for E^i = 0.01, Fig. gives a solution of a 
remnant at 8-10 kpc and an age of 3000 — 6000 yr. The 
association 3C58-SN1181 suggests an explosion energy of 
the order of 10^* erg, roughly ten times lower than the 
lowest inferred. We recall that the association have been 



questioned by Huang (1986), but recently Stephenson fc 



[Green (1999) have pointed out that, on the basis of an 



update of historical information, the association should be 
reliable. If so, the Sedov model does not provide a proper 
description of all the observational evidence. 

In order to understand the interaction between 3C58 
and its environment, and to understand the origin of the 
soft X-ray component, it is useful to compare its radio 
and X-ray emission. The soft X-ray maps of Fig. |^ show 
that the X-ray emission is in any case confined within 2'- 
4' from the center, and therefore the presence of the shell 
is in agreem ent with the lack of detec tion of a radio shell 
at r > 5' of Reynolds fc AUer (1985) . It is interesting to 
note that there is a close correspondence between the ra- 
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dio morphology of the outer regions of the nebula and the 
soft X-ray maps (Fig. ^ upper panel). This is also ob- 



served in other shell-like young SNR like Kepler (Matsui 
et al. 1 384) and Cas A ( [Keohane et al. 1996 ), while, on the 
other hand, the radio emission of the Crab nebula is four 



times greater than its X-ray counterpart. Reynolds & AUer 
(1988)[]have pointed out that the radio image of the rem- 
nant at 1446 MHz shows confined-edge emission at some 
locations of the outer regions of the nebula, and some of 
these also show limb-brightening. According to them, this 
may be explained if the edge of the nebula is sweeping up 
moving material, e.g. ejecta, like in th e "inhomogeneous" 
model of Reynolds fc Chevalier (1984) , or the shock model 
invoked by Sankrit fc Hester (1997) to explain the [O HI] 
emission seen at the boundary of the Crab. This may ex- 
plain the small amount of limb brightening and is also in 
agreement with the relatively slow shock speed measured 
in X-rays. In fact, the expected shock speed in the mov- 
ing ejecta material is ~ 300 km ( Reynolds fc Chevalier 



1984 p|50-20 km s'^ in the Crab according to ^ankrit fc 
Hester 1997), and we observe Vg = 450 km s^^. Therefore 



the thermal component we observe in the outer rim of 
3C58 may be associated with the expansion of the nebula 
in the inner ejecta core. Note that this is also in agree- 
ment with the discrepancy between the large speed mea- 
sured in 3C58 filaments {'^ 900 km s~^, Fesen 198^ ) and 
the X-ray derived speed. In fact, the filaments are com- 
posed of material ejected in the explosion itself, pushed on 
and accelerated by the synchrotron nebula, while the X- 
ray emission is due a shock expanding into moving ejecta. 
However, the low value of the X-ray emitting mass would 
imply that the interaction is only at its beginning. 

Finally, it should be noted that the above conclusions 
rely on the hypothesis of fast electron-ion equipartition. 
Unfortunately, it is not possible to independently measure 



the proton temperature Tp of the 3C58 shock. Bocchino 



et al. (|1999)| have shown that < T„ < 2Tp. for a shock 



region of the old Vela SNR, but [Hughes et al. (2000)| found 
evidence for Tp - A5Te in the 1000 yr old SNR EO 102.2- 
7219, suggesting that non-equipartition may be common 
among young SNRs. If also in case of 3C58 Tp is ~ 45 
times higher than the X-ray derived electron temperature, 
than a Sedov solution with Eq ^ 3 x 10'^^ erg, would be 
compatible with a distance of 3.2 kpc and the association 
with SN1181. 



4. Summary and conclusions 

The high spatial resolution of the EPIC camera on-board 
XMM-Newton has allowed us to look with unprecedented 
detail at the X-ray emission from the fiUed-center SNR 
3C58. We have presented soft (0.1-1 keV), medium (1- 
2 keV) and hard (2-10 keV) energy images, and we have 
found a close correlation of the soft image with radio emis- 
sion at 1446 MHz. A quantitative study of the source pro- 
file shows the effect of synchrotron burn-off and indicates a 
different slope of the soft X-ray radial profile compared to 
that obtained in the hard band. Spatially resolved spectral 



analysis carried on in rings with Ar ~ 8" has allowed us 
to derive the relation between the spectral index and the 
distance from the core. In contrast with previous ASCA 
results, we have not found evidence for thermal black-body 
emission from a central source, and we have placed an up- 
per limit to this component of Lx = 1-8 x lO'^'^ erg s~^. 
This stringent upper-limit rules out thermal emission from 
the whole surface of the putative neutron star at the cen- 
ter of 3C58, and also rules out the "outer-gap" model for 
hot polar caps of Cheng et al. (1986), unless the pulsar 
has very low rotational energy losses. 

The large effective area of the XMM-Newton mirrors 
has allowed us to perform, for the first time, a spectral 
analysis of the outer edge of the 3C58 X-ray nebula, ex- 
cluding the inner and brighter regions. We have found 
that non-thermal emission is still responsible for most of 
the emission there, but a soft thermal component is re- 
quired to better fit the spectrum below 1 keV. If repre- 
sented with an optically thin emission model, this com- 
ponent gives kT=0.2-0.3 keV. If it is associated to the 
Sedov expansion of the 3C58 shell, it is incompatible with 
the association between 3C58 and SN1181, unless there 
is strong deviation from electron-ion equipartition. If it is 
due to the expansion of the nebula into the inner core of 
moving ejecta, the X-ray spectra characteristics and the 
radio morphology of the outer nebula can be more easily 
reconciled. 
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Fig. 2. MOSl+2 images binned at pixels of 2" and 
smoothed using a gaussian with cr = 3 pixel to enhance 
the morphology of the faint nebula edge. Upper panel: soft 
image (0.1-1.0 keV). Center panel: medium energy image 
(1.0-2.0 keV). Lower panel: hard image (2.0-10.0 keV). 
Overlaid are contours representing the weakest (3cr = 0.15 
Jy beam~^) 1446 MHz radio contour and the one corre- 
sponding to a flux density 4 times higher (both in white, 
from Reynolds fc AUer 1988| ), and six X-ray contours (in 
yellow, from 1/100 of the peak to the peak value, logarith- 
mically spaced). 
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Fig. 3. Spatially resolved spectral analysis with a power- 
law model of the X-ray nebula associated with 3C58. We 
show the results obtained on the 8 annuli and the "edge" 
regions (the last point) as defined in the text. Black is 
MOS, red is PN, dotted bars are fit with free Nh, solid 
bars are fit with Nh fixed at 4 x 10^^ cm~^. Regions are 
centered at 2''05'"38.0^ +64'*49'"37" (J2000). A linear fit 
to the MOS points corresponding to fits with Nh fixed at 
4 X 10^^ cm^^ is shown in the lower panel. 
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Fig. 5. Same as in Fig. ^, but here we report the ratio be- 
tween the 0.5-2.0 keV flux due to the additional compo- 
nent (mekal in the upper panel, and bbody in the lower 
panel) to the total 0.5-2.0 keV flux versus the distance 
from the center. Arrows indicate upper limits. Goodness 
of fits are reported in Table |l]. 
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Fig. 6. The plot shows that allowed regions of the B — P 
plane according to the cla ssical model for pulsars (e.g. 
Manchester fc Taylor 1977 ) in solid lines (assuming the 
initial period was much less than the present period) . We 
show the case for an 800 yr old pulsar (i.e. assuming that 
3C58 is related to SN1181) and for a pulsar with an age a 
factor 10 larger. The dashed line represent the B — P val- 
ues allowed by the "outer-gap" polar heated caps model 



of Cheng et al. (1986) and our derived upper limit to the 
black-body luminosity of the central source. The intersec- 
tion of the observational and theoretical curves gives the 
allowed B — P values (thicker part of the isochrone lines). 
If 3C58 is the remnant of SN1181, the model would pro- 
vide an unreasonably high B. Older pulsars and/or low 
rotational energy losses may yield more reasonable B val- 
ues. 



Fig. 7. Loci of allowed SNR age and distance according 
to a simple Sedov analysis of the thermal component of 
outer rim X-ray emission of 3C58. The solid line labeled 
"Geom." gives the loci allowed by the simple geometrical 
relation betwe en the real radius of the shell and the SNR 
age (Eq. 2 of Kassim et al. 199^ ). The radius is linked 



to the distance via the apparent shell radius, and we have 
used 2.5' to derive the relation. The dashed lines represent 
the loci of the solutions allowed by the Sedov relation R = 
UiEsi/noY^Hl^^ where no is derived from the emission 
measure of the mekal component assuming a spherical 
emitting volume with a radius of 120"and a filling factor 
of 25%. The intersection between the dashed lines and the 
solid line gives a solution. Vertical dashed lines mark the 
distance estimate of 2.6 and 3.2, due to Green & Gull 



1982)1 and [Roberts et al. (1993)| , respectively. 



